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Abstract 

Transverse spectra are calculated for various types of hadrons stemming from Au Au collisions at V* - 200 GeV. 
We utilize a quark recombination model based on generalized Boltzmann-Gibbs thermodynamics for local hadron 
production at various break-up scenarios. 
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1. Introduction 

Hydrodynamics turned out to be a successful tool 
in describing the spacetime evolution of the hot dense 
quark-gluon plasma (QGP) created in heavy-ion colli- 
sions (HIC) until rehadronisation Hl-lS]. Various types 
of quark recombination models based on Boltzmann- 
Gibbs (BG) thermodynamics l4ll- lfloll have been used 
with success to reproduce the soft part (pj < 2 GeV) 
of transverse hadron spectra observed in Au Au col- 
lisions at V* = 200 GeV HI -El. However in a 
QGP at hadronization there are long ranged interac- 
tions. Due to these, non-additive effects may not be- 
come negligible in the thermodynamical limit. The ratio 
of non-additive to additive terms in general scales like 
surface X force range /volume. This quantity may not 
vanish in the limit of large volumes. In this case the BG 
thermodynamics has to be gen eralized l30ll . 
Recent RHIC, Tevatron 0] and LHC HH trans- 



calculation. In the stringy M43I1 . and in the distributed 



verse hadron spectra are well fitted by the distribution, 
dN/p T dp T dy = A(l + p T /po)~", at p T < 5 GeV This 
formula has also been successfully applied in statisti- 
cal and financial models. There are various attempts to 
clarify the underlying physics that may result in such 
distributions. In works IU9ll - ll4lll among others super- 
statistics, non-additive entropy or energy formulas due 
to long range correlations or interactions, fractal filling 
patterns in the phase space are discussed. However the 
microdynamical emergence of such a hadron distribu- 
tion in a HIC is still not yet presented in a convincing 
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mass [42] quark coalescence model, quarks are assumed 
to be thermally distributed with the energy distribution 

f(E)=A[l + (q-l)E/Tr 1Kq - 1) . (1) 

Here T is the temperature of the QGP, q is a parame- 
ter showing the deformation of the BG thermodynam- 
ics. Sudden coalescence of n quarks results in hadron 
distributions of the cut power-law form with scaled pa- 
rameters 



n(q H - 1) = {q q - 1). 



(2) 



It is also an open question whether the created mat- 
ter reaches thermalization in the deconfined phase or it 
thermalizes in the hadronic phase, if at all. 
In this paper we show that various hadron distributions 
all point out the hadron energy in the comoving frame as 
scaling variable. We assume that hadrons are produced 
at an instant propertime and the image of this distribu- 
tion in the detector can be calculated via the Cooper- 
Frye (CF) formula. This calculation fits well to all types 
of transverse hadron spectra measured in — 200 
GeV AuAu collisions up to p T = 5 GeV. 

2. Transverse Spectrum Via the Cooper-Frye For- 
mula 

From the local equilibrium hadron distribution the 
transverse spectrum is calculated via the CF formula 



dN 



d 3 N 



dp 3 prdpjdydip 



= j F(u'pi)pi 



<fV. (3) 
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with p' four-momentum of the detected particle, 
local flow four- velocity (normalized to m,m ! = +1) and 
x' spacetime coordinates, parametrized as 



2n 



(rch£, rsh^, rcosa, rsina), 

(ych77, ysh77, yvcos^, yvsin</>), 

(mj chy, ttij shy, pj cos <p, pj sin <p). (4) 



Here t = Vf 2 - z 2 is the longitudinal proper time, r the 
transverse (cylindrical) radius, ( - 1 In '-0 the longitu- 
dinal coordinate rapidity (hyperbolic arc angle) and a 
the asimuthal angle belonging to the spacetime point 
x'. Similarity v is the transverse flow velocity, y = 
1 / Vl - v 2 the associated Lorentz factor, 77 the longitu- 
dinal rapidity of the flow and <p is its asimuthal angle. 
Finally pj is the detected transverse momentum, mj = 



4 



m 2 + p 2 the corresponding transverse mass, while y 
the observed longitudinal rapidity, E = mj cosh(y) the 
energy and <p the asimuthal angle of the detected paticle. 
This way the integration measure in (f3j) is given by 

Pid^cr' = mj cosh(y — f) rr d^drda 

+ mj sinh(y — ( ) r drdrda 

+ pj cos(</j - a) rr drdrda 

+ pj sin(^ - a)T drd^dr. (5) 

Considering a flow, constant and radial in the transverse 
plane (0 = a, v=const.) and scaling in the beam direc- 
tion (77 = £), the energy in the comoving frame is 

E c " = u'pi = 

= y{mj cosh(y - £) - vpj cos(<^> - a)). (6) 

Assuming an instant hadronization at a constant longi- 
tudinal proper-time, 17,, and a continuous evaporation of 
hadrons at the cylindrical surface, the measured spec- 
trum splits to four components: 



d 3 N 

e— = y + y 1 + y 2 + y 4 . 

dp i 



(7) 



These are 



To 



In R +z 

Tb J da J" drr J" d( mj X 

-z 
x cosh(y - £)F(Tb, f , r, a), 

In R T b 

dnn 1 x 
000 

x (sinh(y - z)F(t, +z, r, a) 

- sinh(y + z)F(t, -z, r, a)), 



zn j\ Tb 

= J" da J" drr J" c 



Y 2 



T3 



R 



drrpj X 



-z 

X cos((p - a)F(j, R, a), 

R +Z T b 



dT T pj X 



X (sm((p - 2n)F(T, {, r, 2n) 
- sm(ip)F(T, £ r, 0)). 



(8) 



Here T3 = for F(E co ) is periodic in a, and Tl 1— > 
taking the z h> 00 limit. This way the spectrum can be 
written in a more convenient form, i.e. as a sum of a 
volume and a surface term. After a shift in the integra- 
tion variables £ — > £ + y and « — > a + <p the transverse 
spectrum takes the form 

00 2n 

d 3 N f f 

£— = A J d£ J daF(-/3E co ) x 



x (£ mr cosh f + (1 - £) p r cos a) (9) 



with 



A£ = 



A (1-0 



ft 

In J" rr(r)dr = nR 2 Tt,, 


27T r(r)TdT — jit 2 7 R. 

T(0) 



(10) 



In the first integral of ( TTOb r(r) = tj,, in the second r(r) = 
R. Assuming the BG hadron distribution F(-/3E co ) = 
exp(-/3E co )/Z(J3) we obtain 



, d*N 



A[i; mjK\(fiymj)I§(fiyvp T ) 

(1 -0 PTKo(Pym T )h(J3yvp T )). (11) 



This result is a decreasing nearly exponential function 
of E* = y(mj-v pj) above pj = 2 GeV. However, using 
the hadron distribution (Q]) we utilize its Euler-Gamma 
representation I311- J36ll . Here c = l/(q - 1) and 



dN 



2n 



p T dp T dy y=0 



da X 



-°° 

x mj cosh ( + (1 - ^) pj cos a) 

X exp{-w[l +(q- l)J3E co ]}. (12) 



2 



The integration over f and a leads to the spectrum 



dK 



v=0 



pjdpjdy 



= A 



%m T Gq(pt) + (1 - Opt G 2 (pt) 
(l + (q-l)/3E*y 



(13) 

It differs from the results presented in ll43"ll by the terms 
iiitGq and PtGi. The factors G, show, in fact, a very 
mild pr dependence. 

Go,i(Pt) = V ^-e a >-h K {a y {b ) 

/ - j r(c) 

./ 

(5 - Y)Pym T 



l + (q-l)j3E* J 
(q - \)Pyvp T _ 
1 + (q - !)/?£ 



(14) 



with = Zj/(n - l) 2 L(n + l,z/) 2 , L(n,x) being the La- 
guerre polynomial, and L(n,zj) = 0. 

3. Results 

The first 10 panels in Figure Q] show the ratios of 
meassured and calculated transverse hadron spectra 
dN I pr dpj dy at mid rapidity (y = 0). Boxes denote 
results stemming from BG thermodynamics, triangles, 
dark and light dots denote the results obtained from the 
generalized hadron distribution Triangles belong to 
£ = 0, dark dots to £, = 0.5 and light dots to £ = 1. 
This variable reflects the surface to volume term ratio in 
the hadronization formula (0. It is clear that the appli- 
cability of the BG distribution is limited to a very nar- 
row range, up to 2 GeV, while the distribution ([TJ fits 
better than a factor 1.5 up to 5-6 GeV for all hadron 
types. Furthermore it is likely that most of the particles 
are produced in the instant hadronization process at Tf,. 
The data/theory ratios are closest to one for ^ = 1 which 
is taken in most calculations in the literature. Panel 1 1 
on Figure Q] shows the extracted values for the tempera- 
ture, T, against those of the non-extensivity parameters 
{qH>- In panel 12 on Figure Q] trans verse flow velocities 
(v/c) are plotted vs. the ratios of surface and volume 
yields (£). The cases ^ = 0.5 and £ = 1 are depicted by 
dark and light dots respectively. The fitted temperatures 
and non-extensivity parameters are close for all hadrons 
(T = 51 ± 10 MeV, q = 1.062 ± 7.65 x 10~ 3 ). We 
conclude that the matter created in these HICs is ther- 
malized according to energy shells, although the energy 
distribution is not the BG exponential 11441 14511 . 
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Figure 1: panels 1-10: Ratios ofRHIC 200 GeV AuAu data and calcu- 
lated spectra. Boxes: results of BG thermodynamics; triangles, dark 
and light spots: results of the hadron distribution {T} with f = 0, 0.5, 1 
respectively; Panels 11-12: fit parameters. 
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